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Abstract. We examine the Stark effect (the second-order shifts in the energy 
spectrum due to an external constant force) for two 1-dimensional model quantum 
mechanical systems described by linear potentials, the so-called quantum bouncer 
(defined by V{z) = Fz for z > and V{z) = oo for z < 0) and the symmetric 
linear potential (given by V{z) — F\z\). We show how straightforward use of the 
most obvious properties of the Airy function solutions and simple Taylor expansions 
give closed form results for the Stark shifts in both systems. These exact results are 
then compared to other approximation techniques, such as perturbation theory and 
WKB methods. These expressions add to the small number of closed-form descriptions 
available for the Stark effect in model quantum mechanical systems. 
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1. Introduction 

There are only a handful of model quantum mechanical systems which admit exact 
solutions, and from which information is easily extracted in closed-form expressions. 
The most familiar examples include the Kepler problem (Coulomb potential for the 
hydrogen atom), the harmonic oscillator (for vibrational states), spherical harmonics 
(for rotational motion), and the infinite square well. These examples all allow for 
a relatively complete analysis of the mathematical descriptions of the position-space 
wavef unctions, including closed-form expressions for normalizations, expectation values, 
and even matrix elements for use in applications such as sum rules and perturbation 
theory. 

The so-called 'quantum bouncer' system, defined by the potential 



has been a staple of the pedagogical literature, in both articles [I] - |5] and textbooks [6] - 
[9]. It not only has solutions given by well-studied special functions (Airy functions p^). 
but it can also be analyzed using a variety of approximation methods for comparison 
to exact results. Most importantly, however, this 'academic' problem has received 
renewed attention as the simplest model for the quantum states of neutrons in the 
Earth's gravitational field [llj. The linear potential plus infinite wall has also been 
used to analyze surface-Landau-level resonance data [12] where electrons are trapped 
between the surface potential and a weak magnetic force given hy F = evpH/c. In 
addition, related recent experimental realizations of motion in such a 'bouncer' potential 
also include systems of atoms [I3], Bose-Einstein condensates [H], and even optical 
analogs ('photon bouncing balls') [TS]. These proverbial 'real life' applications help 
to promote the quantum bouncer to the small set of important physical realizations of 
mathematically soluble problems which are accessible to students of quantum mechanics 
at the undergraduate level. 

A simple extension of this problem, the symmetric linear potential, given by 
V{z) = F\z\, also shares Airy function solutions. Moreover, it has the same symmetry 
property as systems such as the harmonic oscillator or symmetric infinite well and so can 
be discussed in a similar context, including using concepts such as parity. For example, 
all three cases can be thought of as belonging to the class of power-law potentials of the 
form 



with A; = 1, 2 and oo for the three systems respectively. One can then systematically 
explore the behavior of energy eigenvalues or even higher order effects, such as 
the Stark shifts [l7j, and how they depend on the power-law exponent (/c), quantum 
number, and physical parameters, such as the mass. 

The Stark effect, the response of the quantized energy levels to the addition of 
a constant external force, is another staple of quantum mechanics texts, with closed 




for < z 
for z < 



(1) 
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form results possible in a number of the same tractable systems mentioned above. For 
example, the addition of a linear potential of the form V{z) = Fz to the harmonic 
oscillator (HO) can be easily accommodated with an exact solution (requiring only a 
simple ^complete the squares^ trick), leading to a second order shift, 

namely a constant, lowering the energies of all of the states. For the symmetric 
infinite well (SIW), with walls at x = ±a, methods involving inhomogeneous differential 
equations [18] (the so-called Dalgarno-Lewis method [19]) and perturbation theory 
summation approaches [20], [Hj have been used to evaluate the second order shift, 
yielding 



1- 15 



(4) 



where 



E^^ ^ ^'"^y X where n = 0,1,2,... (5) 

define the unperturbed energies in an infinite well of width 2a. In contrast to the 
harmonic oscillator, the Stark shift here is negative only for the ground state (as required 
by second-order perturbation theory) but positive for all excited states, which is an 
interesting qualitative difference. It is in this context that we will examine the Stark 
shifts for the quantum bouncer and the related symmetric linear potential. Given the 
simple nature of the Airy function solutions (or rather the differential equation they 
satisfy), it's perhaps not surprising that exact expressions for the Stark shifts in both 
systems can also be derived, and that is the topic of this paper. 

In the next section, we review the solutions of the quantum bouncer problem, 
focusing on providing a list of exact results needed for quantities such as normalizations, 
expectation values, and other matrix elements. We find that the Stark shifts are actually 
trivially obtained in this case. In Sec. [3l we provide the corresponding results for the 
more interesting case of the symmetric linear potential. For this example, we make 
use of the fact that the simple form of the Airy differential equation provides dramatic 
simplifications for the Taylor series expansion of the exact energy eigenvalue condition 
for the 'perturbed' case, and we find closed form expressions for the Stark shifts in this 
system as well, one of our main results. In both cases, we compare the exact results 
to WKB approximations for the Stark shifts, which have been discussed earlier in the 
literature [17], as well as with second-order perturbation theory. Finally we put our 
results into context in Sec. H] where we discuss our conclusions. 
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P=i^Y and /3. = |^ = |^. (8) 



2. Stark effect for the quantum bouncer 

We begin by reviewing the solutions for the quantum bouncer problem, defined by the 
potential in Eqn. ([T]). The Schrodinger equation in the region 2; > reduces to 

- + F^^niz) = EMz) , (6) 

which can be written in the form 

^:{x) = (x - P^)^^{x) (7) 
using the change of variable z = px and the definitions 

\2mF) Fp ~ So 

The solutions of Eqn. ([7]) are the two linearly-independent Airy functions, Ai{x — 
and Bi{x — /?„). (The now standard reference on the application of Airy functions in 
physics is Ref. [lO].) The Bi solution (which diverges for large positive argument) does 
not satisfy the boundary condition ipniz 00) = and so is excluded. The energy 
eigenvalues are then determined by the boundary condition imposed by the infinite wall 
at the origin, namely that iIj{z = 0) = Ai{—f3n)- The quantized energies are then given 
in terms of the zeros of the well-behaved Airy function, Ai{—(n), with En = +Cn^o- 
In many other systems, the explicit manner in which the application of the boundary 
conditions leads to discrete energies is less than mathematically obvious. In this case, 
however, one simply 'slides' the single Ai{x) solution along the axis (via the /?„ shift) 
until one of its zeros coincides with the origin and therefore satisfies the boundary 
condition at the infinite wall. (We note that Eqn. ([6]) can be easily solved in momentum 
space, and then connected to the integral representation of the Airy function |9]. The 
energy eigenvalue conditions for this problem, and especially our systematic expansion 
of them, are much more easily derived directly from the position space form of the 
Schrodinger equation.) 

Standard handbook results for the asymptotic behavior of the Ai{x) zeros [22] give 



En ~ So 



Y(--1/4) 



2/3 



£0 



Svr , , , 

— (2n-l/2) 



2/3 



(9) 



for large n (where the labeling starts with n = 1.) We note that the WKB prediction 
for the energy eigenvalues for this system is 



J 

Jo 



'2m{En- Fz)dz= {n + CL + CR)fnT with n = 0,1,2,... (10) 

/o 

and the appropriate values of the matching coefficients are Cl = 1/2 (at the left turning 
point, where there is an infinite wall boundary condition) and Cr = 1/4 (for the 
smoother potential at the right turning point). Solving for the quantized energies gives 

n 2/3 



EniWKB) 



37r(n + 3/4) 



So (11) 



2 

where n = 0,l,2,..., which agrees with the 'handbook' result in Eqn. ([9]) for large n. 



Stark effect for linear potentials 



5 



The wavefunctions for positive z are then given by 




(12) 



and the orthonormahty properties of these solutions can be investigated very 
straightforwardly. For example, the normalization constant can be determined in closed 
form by using one of the simple identities, Eqn. (15^ . involving integrals Airy functions 
collected in the Appendix, where we find that 

^"^v^^- '''' 

(Gea-Banacloche [23] was the first to find this relation numerically, and soon afterwards 
Vallee |24| and Goodmanson [25] demonstrated it analytically. In this work, we choose 
Ai'{—Cn) in the normalization instead of |74z'(— (^„)|, as done in Refs. [23] - [25], in order 
to simplify some of the expressions below; such a choice of phase has, of course, no 
physical significance.) One can then just as easily confirm that {ipn\'4'm) = for n 7^ m 
by using Eqn. (l58ll . 

As noted by Goodmanson [25], the expectation values of the potential and kinetic 
energies can be evaluated by using the identities in Eqns. (!55l) and (!57j) respectively to 
find 

roo 

Z I'lpniz)]'^ dz 



{n\V{z)\n) 




(14) 



and 



2m 



(n|p^|ra) 




Cn)fdx 



(15) 



where an integration by parts has been used in Eqn. (fT5|) to bring it to the form of the 
integral in Eqn. (|57l) . (One can also make direct use of the Airy differential equation, 
A" = zA, to rewrite {n\p'^\n) in terms of the same integral in Eqn. (IT^ .) These 
calculations confirm that both quantities are consistent with the virial theorem. For 
future reference, we find that the dipole matrix elements are given by 

2p 



{i'n\z\i'k) 



(16) 



(Cn - Ck)' 

where we use Eqn. ( I59l) . 

Turning now to the Stark effect, we note that the result of a perturbing potential 
of the form V{z) = Fz is trivially realized since the form of the solutions remain the 



Stark effect for linear potentials 



6 



same with the global substitution F 
obtained from Eqn. ([8]) as 

2/3 

En = (^oCn 



[F + F). The new energy eigenvalues are thus 



:i7) 



With this simple form, the predictions for the first-, second-, and third-order energy 
shifts are given by a simple expansion, namely 

I (F^ 



eP 



eP = - 




and 




1^ 



The first-order perturbation theory result. 




(19) 



is easily confirmed by the same integral as in Eqn. ( |T^ . The second-order energy term 
is given by the standard expression 

\{n\V{z)\k)\' 



E^' = E 



SO that the predicted Stark shift for the quantum bouncer is 

^(2) 1 




(20) 



(21) 



It has been shown in Ref. [26] that this result (and the corresponding double summation 
giving the third-order expression, E^^) provides new constraints on combinations of 
zeros of Airy functions, which can be easily confirmed numerically using tools such as 
Mathematica® . This is one example where novel applications of standard quantum 
mechanical methods can provide new results in mathematical physics. 

We also note that the WKB approximation in Eqn. (1 10 1) has the same global 
redefinition of F — > F = F + F, yielding the same additional factor of (1 + F/F)^/^ 
found in the exact result in Eqn. flTTl) . so that the WKB approximation gives the correct 
large n behavior of the first- and second-order shifts as well. 



3. The symmetric linear potential 

The case of the symmetric linear potential, a seemingly small variation on the quantum 
bouncer potential, namely 

V{z)=F\z\, (22) 

adds a number of interesting features, most notably making the evaluation of the F- 
dependent shifts in energy levels more challenging. But use of the simplest properties 
of the Airy function solutions leads to a straightforward calculational algorithm which 
can still give closed-form results for the energy shifts to any desired order in F/F. We 
begin by reviewing the solutions of the symmetric linear potential. 
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The energy eigenstates in this symmetric potential can be classified by parity, 
and the odd-parity states are automatically related to those of the 'half-well' problem 
(namely the quantum bouncer above) by 

+M^)/V2 forz>0 
^" ^ ^ ~ I -M-^)/V2 for z < ^^^^ 

where the ipn{z) are given by Eqns. f|T2|) and f|T3|) and the normalization is modified. 
Henceforth, we label the energies for these states as -E.^^^ = Cn^o- 

The corresponding even states must still be of the form -ipl^^^z) = Ai{z/p — 
(at least for 2; > 0), but the appropriate boundary condition is now that ip^z = 0) = 
Ai'{—Pn) = 0, so the energy eigenvalues for these states are given by i?^"*"-' = Xn£o, where 
— Xn are the zeros of the derivative of Ai{x). Using the integral in Eqn. fl5^ . we find 
that the appropriately normalized even solutions can be written in the form 

1 




Xr. 



V2pxn Ai{-Xn 

once again with an analytic result for the normalization. A standard handbook 
result for the zeros of Ai'{x) gives 



(24) 



|(--3/4) 



2/3 



80 



1/2) 



2/3 



(25) 



Comparing this to Eqn. ([9]), we see that the energy eigenstates thus satisfy E^^'* < 

< -E^+i and interleave, as in any Sturm-Liouville problem, and alternate in parity, 

as for any symmetric potential. We note that the WKB prediction for the symmetric 

linear potential is dictated by the quantization constraint 

r+Ew/F I 

' '~ ' (rr+ 1/4+ l/4);i7r with n = 0,1,2,... (26) 



2miE^-F\z\)dz 



since both turning points are now 'soft' [Cl 
"37r(n + l/2)^'/' 



R 



with n = 0,1,2, 
1/4), giving 



So. 



(27) 



2n 



This expression agrees with Eqn. ([9]) (for n = 2n) and Eqn. ( 125|) (for n 
expected. 

Simple calculations then show that virial theorem results in Eqns. (1141) and ([T5 
still hold. Furthermore, the 'diagonal' dipole matrix elements are given by 



as 



because of parity constraints, while using the integral in Eqn. (!59|) we find 



(+)\ 



2p 



(28) 



(29) 



y/XkiXk - Cn)^ ' 

Turning now to the Stark effect, the addition of an external constant field 
V{z) = Fz, changes the potential asymmetrically to 

' +{F + F)z for z > 

-{F-F)z forz<0 ' 



V{z) = V{z) + V{z) 



(30) 
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but the Schrodinger equation still supports Airy function solutions in both regions. The 
most general solution satisfying the boundary conditions at ±00 is now 

I ( \- \ aRAi{+z/pR - En/iFRpn)) for 2; > 
"^"^^^ \ ar.Aii-z/pr.-Ej{Fr^p^)) for z < ^'''^ 

using the definitions 

Fr = F{1 + A) p^ = p(1 + A)-i/3 (32) 

Fl = F{1 - A) pL = p(l - A)-i/3 (33) 

where A = F/F and the R/L notations refer to the z > and 2; < regions respectively. 
The boundary conditions which must now be satisfied at z = are the continuity of 
both ipn{z) and ip'n{z) and applying these we find the eigenvalue condition 

—En \ 



A)2/3 



We note that in the absence of the perturbing field we have A = and the eigenvalue 
condition reduces to 

G(K, 0) = Ai (^) Ai' (^) = (35) 

which is satisfied by either En/So — Xn (giving the even states) or (^^ (giving the odd 
states), as expected. 

Assuming that A << 1, one can expand G{En,A) about the unperturbed values 
using standard Taylor series expansions of the Ai{x) and Ai'{x) functions. This 
expansion is immensely simplified by repeated use of the fact that the Ai{x) satisfies 
the Airy differential equation, namely Ai"{x) = xAi{x). For example, if we need 

Aiiy + 5)^ Aiiy) + 5Ai'{y) + ^^Az"iy) + ^^Ai/"{y) + ^^Ai^''\y) + • • • (36) 

we can use the relations 

Af{y) =yAi{y) (37) 

Af'iv) = [At"iy)]' = [yAtiy)]' = At{y) + yAt'{y) (38) 
At^^^\y) = [A^iy]]' = [Atiy) + yAt'{y)]' = 2Ai'{y) + yAi"{y) 

^2Ai'{y)+y^Ai{y) (39) 



Thus, the derivatives at any order can be expressed as a combination of Ai{y) and 
Ai'{y) (each times polynomials in y), but no higher derivatives. Then, since we are 
expanding about y — —Xn or y — — one or the other of these terms will vanish 
(cither Ai{—(n) = or Ai'{—Xn) = 0). These simplifications leave the expansion of 
G{En, A) proportional to the non-vanishing term squared, multiplying a polynomial in 
A with coefficients containing polynomials in Xn or 
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= xSo{l + RiA + R2A^ + R^A^ + ■ ■ ■) (40) 

where x = Xn or (n, expand G{En, A) in terms of powers of A, and make repeated use of 
the identities in Eqns. ( 1371) - ( HOl) . The resulting expression in powers of A must vanish 
identically, so that its coefficients (which will contain the dimensionless Ri factors) must 
vanish term-by-term. 

As an example, expanding around En = (n^o for the odd states (for which 
Ai{—(n) = 0) gives the constraint 

G{Et\A) = [Ar{-Q)r |-2Cni?iA - (ii±^M^A^ 

i3Rs + 2R^-2RiOCn ^s , ] _o (4^^ 



up to third order in A = F/F. We immediately see that i?i = i?3 = as dictated by the 
symmetry of the problem and perturbation theory results. Recall that the first-order 
shift in Eqn. ( fT9l) vanishes in this case due to parity (since the diagonal dipole matrix 
elements vanish, as in Eqn. fl28l) ). The third-order shift, given by 

;.(3) _ V V {^\nm{k\m\j){j\m\n) _ ^^y^.,. ^ (ao) 

-2^2^ (0) p(0)wp(0) p(0)x 2^ . (0) (0)x2' ^^^> 

can then be seen to explicitly vanish as well for similar reasons. 

The new non-trivial result is that R2 = —7/9, so that the Stark shift for the odd 
states is given by 

4-..._^(|)^C„£„,^-I(|)^<-. (43, 

In exactly the same manner, expanding about En = Xn^Oi we find the corresponding 
result for the Stark shift for the even states, namely 

If we include the perturbing potential, the WKB prediction for the energy 
eigenvalues is now determined by the quantization condition 

n ^En - {F\z\\Fz) dz = (n + l/2)7ih (45) 
where n = 0, 1, 2, ... and the classical turning points are 

a± = ±— ^ . (46) 
The resulting approximation for the energy eigenvalues is found to be 

F„ = F„ I 1 - — ^En[l--A' + ---j (47) 
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where the are given by the F = WKB resuh in Eqn. fl27p . In this approximate 
approach, the first- and third-order terms (in A) do indeed vanish, while the second- 
order Stark shift is given by 

2 



3 



9 I F j " 



(48) 



with no dependence on parity, since that concept plays no role in the semi-classical 
WKB approach. (The Stark shifts in various power-law potentials have been examined 
in Ref. [I?] using just such a WKB approach, where this last result was first noted.) 

Interestingly then, we find that the exact second-order results in Eqns. (j^3|) and 
(jH]) bracket the approximate WKB result and in fact give it as the 'average' effect, but 
they do depend on the parity of the unperturbed eigenstate. Both approaches give the 
correct dependence as being proportional to the unperturbed energy eigenvalue. 

Finally, in the spirit of Ref. [26] , we can use the standard second-order perturbation 



theory expression in Eqn. 
to find 



20|) and the dipole matrix elements in Eqns. (1281) and fl29l) 



-4 



E 



Xk{.Xk - CnY 



4 Sn 



XniCk Xn) 



which leads to two new constraints on the Xn, Cn, namely 




5Xn 

36 



and 



E 



7Cn 

36 



(49) 
(50) 

(51) 



iCk - XnY 36 Y Xkixk - CnY 

both of which are easily confirmed numerically using programs such as Mathematica® . 
(In fact, in the spirit of Ref. [23], we actually found numerical evidence for these 
expressions first, and then derived them in closed form using the approach followed 
above.) 

We note that the expansion of G{En, A) can be continued to essentially arbitrarily 
high order, by extending the Airy function identities in Eqns. (137|) - (HOl) . Comparing 
the next non-trivial result {E^'^^) to standard expressions for the 4th-order energy shift 
(see, e.g., Ref. [27] for the form), for example, yields a complicated constraint on the 
(n, Xn in the form of a triple infinite summation. 



4. Discussion and conclusions 



We have derived closed form expressions for the second-order energy shift due to a 
constant external field (Stark effect) for two model systems, using simple properties 
of the solutions of the Schrodinger equation for the linear potential relevant for both 
problems. These expressions add to the handful (almost literally) of exact results for the 
Stark shifts which are available for mathematically tractable model systems in quantum 
mechanics (including the hydrogen atom, harmonic oscillator, and infinite well.) 
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The results in Eqn. fl43l) and fl44l) for the symmetric hnear potential are very similar, 
but not identical to, a simple WKB prediction, but with an explicit dependence on the 
parity of the state. In contrast, the exact results for the quantum bouncer are trivially 
implemented and agree with both the WKB approach as well as perturbation theory 
(up to at least 3rd order [26].) 

We can now compare the results for the symmetric linear potential to two other 
familiar cases, such as the harmonic oscillator result in Eqn. ([3]) or the infinite well in 
Eqn. (jll) where the first order results both vanish as well. The form of the second-order 
shift in Eqn. ( 1201) requires that the result for E^^ for the ground state of any system 
(assumed non-degenerate) is always negative and that is indeed confirmed. Similarly 
to the harmonic oscillator, the second-order shifts for all higher states here are also 
negative, but are proportional to the initial energies, and not constant. This result is 
conceptually consistent with the WKB discussions in Ref. [17] and has the dependence 
on quantum number predicted there, but the simple WKB approach does not correctly 
distinguish between the results for different parity quantum solutions. 

One well-known textbook author [28] in discussing second-order perturbation 
theory has noted that ...the response of any system to a perturbing force is to deform at 
the expense of that force and hence to decrease the potential energy of interaction^^ . He 
concludes, however, that ^\..on physical grounds we expect the second-order correction 
to lower the energy in general but have proved this to be the case only for the ground 
state." Our analysis provides other examples of where this intuition is proved correct, 
with the explicit form derived in closed form for all quantum states in the system. 

Motivated by such observations, the discussion in Ref. [T7] used WKB 
approximations to examine the 'cross-over' behavior of power-law potentials of the form 
in Eqn. ([2]) where for small k {k = 1,2 for the symmetric linear potential and oscillator) 
the second-order shifts were all negative, while for large /c (/c — *• cxd for the infinite well) 
all values of E^^ were positive, except as required for the ground state. One wanted to 
see just how 'stiff' the potential had to be for the large n (or even n > 1) behavior of 
the second-order shifts to change qualitatively in sign. The results derived here provide 
additional evidence of the validity of such a WKB approach to that general question, 
but also give more examples of exact closed-form expressions for the Stark effect, and 
how they depend on the details of the quantum solutions, in this case their parity. 

In our review of the properties of the quantum bouncer and symmetric linear 
potential, we have collected several existing closed-form results related to the two 
systems, while providing a number of new analytic examples. Each of these results 
depend on relatively straightforward identities (such as the integrals in the Appendix) 
or other properties of the Airy functions, all directly related to the simple form of the 
differential equation they solve. Many of the 'standard' model quantum mechanical 
systems which are of relevance to important experiments (Coulomb potential, harmonic 
oscillator, spherical harmonics) require moderately sophisticated mathematical physics 
techniques to obtain even some of the simplest results; one exception is the infinite 
well, which relies only on simple integrals involving familiar trigonometric functions. 
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Another often discussed problem is the Stark effect for the plane rotator [22] - for 
which the application of perturbation theory is less than straightforward, so much so 
that a number of textbooks have evidently gotten it wrong [33] - [37]. The correct 
results in that case can also be obtained by explicit expansion of the exact eigenvalue 
condition, but this requires knowing the properties of Mathieu function. In contrast 
to these cases, the results shown here are all derivable (or can be easily checked by 
simple differentiation) using only the fact that Ai"{x) = xAi{x) and little else, making 
them very accessible to students of quantum mechanics at the undergraduate level, 
without the need for advanced mathematical background. The fact that they satisfy 
the Schrodinger equation is all that's needed. 

Our focus on these analytic solutions is motivated by recent discussions of new 
mathematical results [23] - [25] involving Airy functions using just such simple identities. 
Many of these results have been in the research literature for decades [38], [39], but are 
only now finding use because of the application of this once 'abstract' problem to novel 
physical systems [H] - [T5] . 

As recently as 1992, one author [5j was led to describe the analytic solutions of the 
quantum bouncer in the following way: ''This technique, which leads to Airy functions, 
tends to be algebraically tedious and not particularly instructive; it is best reserved for 
more important potentials, such as the harmonic oscillator or Coulomb potentials." 
Given the increased physical relevance of this model system, and its mathematically 
tractable nature, we argue both for its utility as a useful pedagogical example and 
for an enhanced position in the pantheon of canonical quantum mechanical showcase 
problems. 

5. Appendix 

We collect here some of the useful identities involving integrals over Airy functions 
needed in this paper. Versions of these results have been derived in Refs. [IE] and [39] . 
In each case, the identities are most easily proved by direct differentiation of the right 
hand sides. (For example, Albright [39] notes that ''Like most results in integral calculus, 
they are easier to check than to obtain") 

We first assume that A{x) and B{x) are any two linearly-independent solutions of 
the Airy differential equation, 

A{x - P) = aAi{x - p) + bBi{x - p) and B{x - P) = cAi{x - p) + dBi{x - p) (52) 
where 

A"{x) = {x- P)A{x - p) and B"{x - P) = {x - p)B{x - P) (53) 

where we consider the shifted arguments needed for solutions of the quantum bouncer 
and symmetric linear potential. We then have the following identities involving indefinite 
integrals of products of the A, 5, their derivatives, and various moments. 
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The first identity below is needed for tlie normalization of wavefunctions, while the 
second and third are used to evaluate expectation values. 

A{x - f3) B{x -(3)dx={x- (3)A{x - (3)B{x -(3)- A'{x - (3)B'{x - f3) (54) 



X A{x - (3) B{x -P)dx = ^(x^ + px- 2p^)A{x - P)B{x - 13) 



+ - {Aix - (3)Bix -(3) + Aix - (3)B'ix - (3)] 
6 

-\{x + 2P)A'{x- P)B\x- (55) 



/ Aix - f3) Bix -f3)dx = — (Sx^ + [3x^ + Af3'^x - 8f3^ - 3)A(x - (3)B(x - (3) 
J 15 

+ — (3x + 2(3) {A'{x - f3)B{x -(3) + A{x - (3)B'{x - (3)} 
15 

- —(3x2 ^ ^ 8/3^)A'{x - /3)B'{x - (3) . (56) 
15 

The following identity involving derivatives is useful for evaluation of the 
expectation value of the kinetic energy operator and can be derived from one in Ref. [39] 
by a simple change of variables. 

J A\x - p) B\x -(3)dx=^ [{A'{x - (3)B{x -(3) + A{x - (3)B'{x - (3)] 

+(x - P)A'{x - f3)B'{x - /5) - (x - p)^A{x - [3)B{x - (3)] . (57) 

We often need to consider integrals involving solutions with two different 'shifts' 
(/3i 7^ (32)1 for example in the confirmation of the orthogonality of eigenstates (first 
identity) or the evaluation of dipole matrix elements, as in Eqns. ( 1161) and fl2^ . 



j A{x - (3i) B{x - (32) dx = ^ {A'{x - (3,)B{x - P2 



{(32 -(3i) 

-A{x-P,)B'{x-P2)} (5^ 



X A{x - A) B{x - P2) dx = i^lp^l-^A{x - [3,)B{x - (32) 

[Pi - (h)^ 



+ 775 ^ {^'(^ - - P2) - Aix - (3i)B\x - (32)} 



(/52-A) (/52-/3i)3 
2 



,A'(x - (3i)B\x - (32) . (59) 



(A - P2)- 
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